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Abstract 



Using a high-statistics lattice simulation of the Ising limit of (A$ 4 )4 theory, we have measured 
the susceptibility and propagator in the broken phase. We confirm our earlier finding of a 
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discrepancy between the field re-scaling implied by the propagator data and that implied 
by the susceptibility. The discrepancy becomes worse as one goes closer to the continuum 
limit; thus, it cannot be explained by residual perturbative effects. The data are consistent 
with an unconventional description of symmetry breaking and "triviality" in which the re- 
scaling factor for the finite-momentum fluctuations tends to unity, but the re-scaling factor 
for the condensate becomes larger and larger as one approaches the continuum limit. In the 
Standard Model this changes the interpretation of the Fermi- const ant scale and its relation 
to the Higgs mass. 



1 Introduction 

The well-established fact that the continuum limit of (A$ 4 )4 theory is "trivial" |l]-0 is 
traditionally interpreted using renormalized perturbation theory. In that picture, the ratio 
M%/v R of the physical Higgs mass to the physical vacuum value vr is proportional to the 
renormalized coupling Xr and tends to zero as l/ln(cutoff). Since vr is fixed by the Fermi 
constant v R ~ 1/(GfV%) ~ (246GeV) 2 , the physical Higgs mass is driven to zero in the 
continuum limit. To avoid that phenomenologically disastrous outcome one must invoke a 
finite cutoff A. A smaller cutoff value allows a larger Higgs mass, but the requirement that 
A > 2M h implies an upper bound on M h ||. 

In general, Vr differs from the bare "Higgs condensate" value measured on the lattice, 
vb = ($)■ The two are related by a re-scaling factor Z: 

- Vb n\ 

VR = 7f (1) 

In the conventional picture Z = Z prop where Z prop is the wavefunction-renormalization 
constant for the propagator of the shifted field $(x) — (<&). Due to 'triviality' one ex- 
pects Z pmp — > 1 in the continuum limit, in agreement with the perturbative prediction 
Zprop = 1 + C(Ar). Lattice data for the shifted-field propagator confirm that Z pTop is quite 
close to 1 0. 

However, a different interpretation of "triviality" has been proposed |9|-3 1| in which there 



are two distinct "Z's" in the broken-symmetry phase. The field $(x) must be divided into 
a finite-momentum piece and a zero-momentum (spacetime-constant) piece, ip. The former 
re-scales by the usual wavefunction-renormalization factor Z prop , but the latter re-scales by 
a different factor, Z v . It is Z = Z v that is needed in Eq. ([!]). Z v is determined by requiring 
the physical mass to match the second derivative of the effective potential with respect to 
the renormalized p>R. 
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where \ is the zero-momentum susceptibility. 



Refs. P-pH argue that in a "trivial" theory the effective potential should be effectively 



given by the sum of the classical potential and the zero-point energy of the shifted fluctuation 
field, which behaves as a free field. This leads to a Veg that is extremely flat in terms of 
the bare field, implying a logarithmically divergent Z v of order ln(A/M/i). Therefore this 
interpretation of "triviality" predicts that Z prop — > 1 and Z^ — > oo in the infinite-cutoff limit. 



A direct test of the "two Z n picture was reported in our previous paper |T2 . There we 
found a discrepancy between the Z prop obtained from the propagator and the Z^ obtained 
from MfrX- Absolutely no sign of such a discrepancy was found in the symmetric phase |I2], 



as expected. Here we report a substantially refined calculation; it involves larger lattices and 
a tenfold increase in statistics. Our previous result is confirmed. 



2 The lattice simulation 



The one-component (A$ 4 )4 theory 



*=£ 



\ £(<&(* + e M ) - $(x)) 2 + ^$ 2 (x) + ^$ 4 (x) 
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becomes in the Ising limit 

Slsing = ~K ^ ^2 ^( X + ^v)^) + ^( X ~ ^)^( X )] 
x )i 

with $(x) = \f2K,(f}{x) and where <j){x) takes only the values +1 or —1. 

We performed Monte-Carlo simulations of this Ising action using the Swendsen-Wang 
I3f cluster algorithm. Statistical errors can be estimated through a direct evaluation of 



the integrated autocorrelation time |[L4|| , or by using the "blocking" [|j|] or the "grouped 



jackknife" [jToj] algorithms. We have checked that applying these three different methods we 
get consistent results. 

We have measured the following lattice observables: 

(i) the bare magnetization, vb = (\<&\), where $ = Y2 X ^( X )/L A is the average field for 
each lattice configuration, 

(ii) the zero-momentum susceptibility 

x = L*[(m*)-(mf], (5) 

(iii) the shifted-field propagator 
G (P) = (j2exp(ipz)(®(x)-VB)($(0)-v B )\ , (6) 

where p^ = ^n^ with n^ being a vector with integer-valued components, not all zero. We 
shall compare the lattice data for G(p) to the 2-parameter formula 

GUP) = -J^-, (7) 



P 2 + m Lt ' 

where mi a tt is the mass in lattice units and p^ = 2 sin ^. If "triviality" is true, then this form 
should give a better and better description of the lattice data as we approach the continuum 
limit; also, Z pmp should tend to unity. 

Another way to determine the mass is to use the method of "time-slice" variables de- 
scribed in ref-lflTll (see also |2l| pp. 56). To this end let us consider a lattice with 3-dimension 
L 3 and temporal dimension L t and the two-point correlator 

d(t, 0; k) = (S c (t; k)S c (0; k) + S s (t; k)S s (0; k)) conn , (8) 



where 

S c (t;k) = ^X>(x,t)cos(k.x), (9) 



^ s (t;k) = -i^0(x,t)sin(k-x). (10) 

X 

Here, t is the Euclidean time; x is the spatial part of the site 4-vector x M ; k is the lattice 
momentum k = (2rr/L)(n x ,n y ,n z ), with (n x ,n y ,n z ) non-negative integers; and (...) conn de- 
notes the connected expectation value with respect to the lattice action, Eq. (f|). In this 
way, parameterizing the correlator C\ in terms of the energy Uk as 

d(t, 0; k) = A [exp(-Wfct) + exp(-c; fc (L t -t))], (11) 

the mass can be determined through the lattice dispersion relation 

3 

mT S (k) = 2(coshcj fc -l) -2^ (l-cosife M ). (12) 

In a free-field theory mTs is independent of k and coincides with mi at t from Eq. (^). 



3 Numerical results: symmetric phase 

As a check of our simulations we started our analysis at K — 0.0740 in the symmetric phase, 
where high-statistics results by Montvay and Weisz [17 1 are available. In Fig. 1 we report 



the data for the scalar propagator suitably re-scaled in order to show the very good quality 
of the fit to Eq. (0). The 2-parameter fit gives m latt = 0.2141(28) and Z prop = 0.9682(23). 
The value at zero-momentum is defined as Z v = mf^x = 0.9702(91). Notice the perfect 
agreement between Z v and Z prop . 



In Fig. 2 we show the values of the time-slice mass Eq. ([12]) at several values of the 
3- momentum and the corresponding result of Ref. [[H]]. The shaded area corresponds to the 
value mi att = 0.2141(28) obtained from the fit to the propagator data. We see that m TS is 
indeed independent of k and agrees well with m\ att . 

Thus, our analysis of the symmetric phase is in good agreement with Ref. [T7| and 
shows the expected "trivial" behaviour. Note that our result for Z prop ~ Z v is in excellent 
agreement with the 1-loop renormalization group prediction Z pert = 0.97(1) f|. 



4 Numerical results: broken phase 

We now choose for k three successive values, k = 0.076, 0.07512, 0.07504, lying just above the 



critical k c ~ 0.0748 [17|. Thus, we are in the broken phase and approaching the continuum 



limit where the correlation length £ becomes much larger than the lattice spacing. 

We used lattice sizes L 4 with L = 20, 32, 32, respectively. This should ensure that 
finite-size effects are sufficiently under control, since in all cases L/£ > 5 fTT|,|T8|j. Finite- 



volume tunneling effects IJJl^Ol should also be negligible, as we explain in the appendix. 
Finally, we also repeated the measurements for k = 0.076 on an L = 32 lattice (with less 
statistics) to confirm directly the absence of finite-size effects in this case. After discarding 
10K sweeps for thermalization, we have performed 500K sweeps (ten times more than in our 
earlier calculation [12|]); the observables were measured every 5 sweeps. Our results for the 
magnetization and the susceptibility are reported in Table 1. We note that at k = 0.076 our 



values are in excellent agreement with the corresponding results of Jansen et al |T9| 



The data for the re-scaled propagator are reported in Figs. 3 — 5. Unlike Fig. 1, the fit 
to Eq. (0), though excellent at higher momenta, does not reproduce the lattice data down 
to zero-momentum. Therefore, in the broken phase, a meaningful determination of Z pTop 
and mi a tt requires excluding the lowest-momentum points from the fit. Fig. 6 shows how 
the chi-squared per d.o.f. of the fit improves, and the fitted value of mi at t stabilizes, as low- 
momentum points are excluded from the fit. Our numbers for mi att and Z prop are reported 
in Table 2. Z prop is indicated in Figs. 3 — 5 by a dashed line. 

The fitted Z pmp is slightly less than one. This fact is attributable to residual interactions 
since we are not exactly at the continuum limit, so that the theory is not yet completely 
"trivial." This explanation is reasonable since we see a tendency for Z prop to approach unity 
as we get closer to the continuum limit. Moreover, we find good agreement between our re- 
sult, Zp r0 p = 0.9321(44), and the Liischer-Weisz perturbative prediction Z pcrt = 0.929(14) || 
at K = 0.0760. The comparison Z pvop = 0.9566(13) with Z pcrt = 0.940(12) at « = 0.07504 is 
also fairly good. 

The quantity Z v is obtained from the product mf att % and is shown in Figs. 3 — 5, as a 
point at p — 0. According to conventional ideas Z^ should be the same as the wavefunction- 
renormalization constant, Z prop , but clearly it is significantly larger. Note that there was no 
such discrepancy in Fig. 1 for the symmetric phase. 

Our data show that the discrepancy gets worse as we approach the critical k. Fig. 7 
shows that Z^ grows rapidly as one approaches the continuum limit (where m\ &tt — > 0). 
Thus, the effect cannot be explained by residual perturbative O(Ar) effects that might cause 
G(p) to deviate from the form in ([I]); such effects die out in the continuum limit, according 
to "triviality." 

The results accord well with the "two Z 11 picture in which, as we approach the continuum 
limit, we expect to see the zero-momentum point, Z^ = mj^x, become higher and higher. In 



the continuum limit the effect should be confined to the zero-momentum point, but it seems 
that, away from the continuum limit, the effect "spills over" into the very low momentum 
modes. This "spillover" would explain why the propagator deviates from free-field form at 
low p, necessitating the exclusion of the lowest few p points to get a good fit to (|7]). As 
we approach the continuum theory the deviation from the Z prop /{p 2 + m 2 att ) form should 
become concentrated in a smaller and smaller range of p, but with a larger and larger spike 
at p = 0. In this limit the shape should become a discontinuous function that is infinite 
at p = and equal to 1 for all non-zero p. The sequence of pictures in Figs. 3 — 5 is quite 
consistent with this expectation. 

The time-slice mass mxs(k) also shows distinctive behaviour at low momentum, as seen 
in Fig. 8 for k = 0.076. Except at low momentum the time-slice mass agrees well with mi a tt- 
However, the time-slice mass at zero momentum, mTs(O), is very significantly smaller than 
iTiiatt (see Table 3). At K = 0.076 our result for mTs(0) is in very good agreement with the 
corresponding result in Ref. |H|. The ratio m - is shown in Fig. 9 as a function of mi att . 
The trend is inverse to that in Fig. 7 for Z v . 

Our interpretation of the data regards m\ at t as the true particle mass and views mTs(0) as 
a symptom of the distinct dynamics of the zero-momentum mode. An alternative interpreta- 
tion might be to regard mxs(0) as the true mass. However, we see two serious objections to 
that interpretation. Firstly, it would entail another "Z" factor, Zq = xmxs(O) 2 , which turns 
out to be roughly 0.88 in all three cases. This Zq, unlike our Z pTop , would not agree well with 
the Liischer and Weisz predictions. (In fact such a discrepancy in the k = 0.076 case was 



noted by Jansen et al |19"|.) Secondly, if mTs(k) depends sensitively on k, as it does near 



k = 0, then the dispersion relation is not well approximated by the (lattice version of the) 
usual form ujk = Vk 2 + const.; in that case the very notion of "mass" becomes problematic. 
At larger k, where mxs(k) does become insensitive to k, it agrees well with mi a tt- 



5 Conclusions 



Our data clearly reveal unconventional behaviour at and near zero momentum in the broken 
phase. Well away from zero momentum the propagator data fit the "trivial" form Z prop /(p 2 + 
m 2 att ), with Z prop in accord with perturbative predictions ||. However, the data deviate from 
this form at very low momentum, presaging an even more dramatic difference between Z prop 
and Z v = ?n 2 att x, where \ is the zero-momentum susceptibility. This difference gets rapidly 
larger as one gets closer to the continuum limit, mi a tt —* 0; see Fig. 7. Because of this fact 
it is not possible to explain the effect in terms of residual perturbative interactions (even 
if one were to use a different value of the mass); perturbative effects should die out as one 
approaches the continuum limit, according to "triviality." Finite- volume effects, including 
tunneling, should be negligible (see appendix). We note that at k = 0.076, L = 20 our 



measurements are in complete agreement with Jansen et al |19| wherever a direct comparison 
is possible. 



Our results are consistent with the unconventional picture of "triviality" and spontaneous 
symmetry breaking of refs. [PHTT|j in which Z^ diverges logarithmically, while Z prop — > 1 in 



the continuum limit. In this picture the Higgs mass M^ can remain finite in units of the 
Fermi-constant scale Vr, even though the ratio M^/vb goes to zero. The Higgs mass is then 
a genuine collective effect and M\ is not proportional to the renormalized self-interaction 
strength. If so, then the whole subject of Higgs mass limits is affected. In view of the 
importance of the issue, both for theory and phenomenology, we hope and expect that our 
lattice results will be checked (and/or challenged) by other groups. 
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Appendix: Tunneling effects 

In this appendix we argue that tunneling effects are negligible in our data. Finite-volume 



tunneling effects in the Ising model have been studied very thoroughly by Jansen et al |19|,|2[ 
In finite volume there is mixing between the +v and — v 'ground states' and the true ground 
state is the symmetric combination with energy Eq s = 0. The antisymmetric combination is 
slightly higher in energy by an amount E^ a which depends on the lattice size L: 

Eo^CL^expi-aL 3 ). (13) 

This formula comes from a semiclassical instanton-type calculation, which also yields theo- 
retical formulas for C and a (Ref. |T9|], Sect. 4.3). For the case k = 0.076, Jansen et al also 
measured E 0a on small lattices (L < 10) P0[ and obtained values a = 0.00358(2), C = 0.101 
that agree well with their theoretical formulas. Using their results we can estimate Eo a 
in our three cases (« = 0.076,0.07512,0.07504 with L = 20,32,32, respectively). We find 
E 0a = 2 x 10~ 13 ,3 x 10~ 8 , and 10" 5 , respectively. 

For the first excited state (in the zero 3-momentum sector) there are also two nearly 
degenerate states, E ls and E\ a . The difference E la — E ls is considerably larger than E 0a . 
Ref. [p(|| 's data for k = 0.076 shows that the ratio of E Xa — E ls to E 0a is around 6 for 
L = 8 and around 4 for L = 10. It therefore seems reasonable to assume that, for larger 
lattices, E\ a — E\ s is not more than 4 times Eq u . (For k = 0.076, L = 20 this implies 
E\ a — Ei s ~ 10~ 12 , so we are fully justified in averaging Ref. fill's results E\ s = 0.3914(12) 
and Ei a = 0.3909(14) in order to compare with our m,Ts(0) result in Table 3.) 

Tunneling effects lead to a modification of the formulas for x and the two-point correlator 



G\. However, the effect of these modifications is tiny in our cases. When expanded for small 
E 0a , the non-trivial factor involved in x ( see Ref. ||19|| ; Eq. (24)) is seen to differ from 
unity only by an amount of order E^JT 2 , where T is the lattice's temporal size. (There is a 
not-immediately-obvious cancellation of the terms linear in E 0a .) This difference from unity 
is tiny (10~ 23 , 1CT 12 , 10~ 7 , in our three cases). Therefore we believe that tunneling effects 
are too small to show up in our data. 
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Ref. L sizc # sweeps « ^ 



2k 



Our data 20 4 5 x 10 5 0.076 0.3015(1) 37.71(22) 
Ref. 0] 20 4 7.5 x 10 6 0.076 0.30158(2) 37.85(6) 



Our data 32 4 2.5 x 10 5 0.076 0.3015(1) 37.70(31) 
Our data 32 4 4 x 10 5 0.07512 0.1617(1) 193.1±1.7 
Our data 32 4 5 x 10 5 0.07504 0.13822(12) 293.38±2.86 
Table 1: The lattice data for the magnetization and the susceptibility. 



L size # sweeps k m latt Z piop Z v 

20 4 5 x 10 5 0.076 0.42865(456) 0.9321(44) 1.0531(232) 

32 4 2.5 x 10 5 0.076 0.42836(500) 0.9312(27) 1.0516(260) 

32 4 4xl0 5 0.07512 0.20623(409) 0.9551(21) 1.2340(502) 

32 4 5 x 10 5 0.07504 0.17229(336) 0.9566(13) 1.307(52) 

Table 2: Our values for mi att and Z prop as obtained from a fit to the lattice data for the 
propagator (see Figs. 3 — 5). We also report Z^ = mf att x where the x's are given in Table 1. 



Ref. L siz0 # sweeps k m TS (0) m latt 

Our data 20 4 5 x 10 5 0.076 0.388(10) 0.42865(456) 
Ref. |H| 20 4 7.5 x 10 6 0.076 0.3912(12) 
Our data 32 4 4 x 10 5 0.07512 0.1737(24) 0.20623(409) 
Our data 32 4 5 x 10 5 0.07504 0.1419(17) 0.17229(336) 
Table 3: The lattice data for the zero-momentum time-slice mass compared with the mass 



from the finite-momentum propagator. The result of Ref. [19|] is an average of the two values 
Ei s an d Ei a reported in their Table 4. 
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FIGURE CAPTIONS 

Figure 1. The lattice data for the re-scaled propagator at K — 0.0740 in the symmetric phase. 
The zero-momentum full point is defined as Z v = m^x- The dashed line indicates 
the value of Z prop . 

Figure 2. The data for the time-slice mass Eq.(12) at different values of the 3-momentum. The 
shaded area represents the value mi att = 0.2141(28) obtained from the fit of the prop- 
agator data. The black square is the result of Ref. [|i~7l . 



Figure 3. The lattice data for the re-scaled propagator at k = 0.076. The zero-momentum full 
point is defined as Z v = mf att x- The very low momentum region is shown in the inset. 
The dashed line indicates the value of Z prop . 

Figure 4. The same as in Fig. 3 at k = 0.07512. 

Figure 5. The same as in Fig. 3 at k = 0.07504. 

Figure 6. The values of the reduced chi-square and of the fitted m\ at t as a function of the number 
of low-momentum points excluded from the fit to Eq.(^). 

Figure 7. The two re-scaling factors Z prop and Z v as a function of mi att . (The continuum limit 
corresponds to mi att —> 0.) 

Figure 8. The time-slice mass mTs(k) for several values of the spatial momentum. The black 



square at zero momentum is the result of Ref. |19| . The shaded area represents the 



value mi a tt — 0.42865(456) obtained from the fit to the propagator data. 
Figure 9. The ratio m ' as a function of mi a tt- 

m latt 
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